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A theoretical study of the dynamical dielectric response of bulk lead at low energies is presented. 
The calculations are performed with full inclusion of the electron band structure calculated by means 
of a first-principles pseudopotential approach. The effect of including the spin-orbit splittings in 
the band structure is analyzed, together with dynamical exchange-correlation and local-field effects. 
Excitations are studied in detail, showing the presence at small momenta of acoustic-like dispersing 
excitations. The character of these modes with sound-like dispersion depends on the direction of the 
momentum transfer q. In general results show strong anisotropy effects on the dielectric response 
of bulk Pb. Comparison with available optical experimental data shows good agreement. 



PACS numbers: 71. 20. Be, 71.45. Gm 

I. INTRODUCTION 

One of the characteristic of lead is that it presents the 
second-highest critical temperature of all bulk elemental 
superconductors (T c = 7.23 K), which has been shown to 
be related to its large electron-phonon coupling (EPC) 
constant. In addition, spin-orbit (SO) coupling has a 
strong impact on the electron-phonon interaction in bulk 
Pb, increasing its strength as much as 44%. 1 Studying 
the electronic collective excitations near the Fermi level 
(Ep) and the SO effects on them would complete the 
description of low-energy dynamics in bulk lead, and it 
is done in the present work. 

For many years the electron-density response of solids 
was studied using a free-electron gas (FEG) model in 
which the electron valence density is parametrized by 
a single quantity: the density parameter r 8 , which 
stands for the average inter-electron distanceP The FEG 
model gave insight on basic properties of the momentum- 
and frequency-dependent dynamical dielectric response. 
However, band structure effects that are missing in a 
FEG model frequently produce dramatic impact on the 
dynamical dieletric response of solids. In particular, in- 
terband transitions (not present in a FEG model) give 
rise, for instance, to a strong red shift of the Ag plas- 
mon frequenc^ or to a negative momentum dispersion of 
the plasmon in CiPl Additionally, interband peaks domi- 
nate the energy-loss landscape in the low-energy transfer 
regime. Also, it has been shown that band structure ef- 
fects play an important role in the screening of charges 
at noble metal surfaces!^ 

For three-dimensional (3D) solids the FEG model pre- 
dicts the existence of a r s -dependent threshold for col- 



lective excitations. Thus, according to the FEG theory, 
plasmons can not participate in the low-energy dynami- 
cal processes near the Fermi surface. However, PineJ^in 
the fifties predicted the existence of a very-low energy ex- 
citation, which should be present in systems with several 
energy bands crossing the Fermi level with different Fermi 
velocities (vf), as it is the case of bulk Pb. This very-low 
energy mode presents an acoustic-like dispersion at small 
momenta: uoap = v ■ q, where v is the constant group ve- 
locity of the acoustic plasmon. Thus, uoap tends to zero 
as q — >• 0. Exchange of acoustic plasmons have been 
suggested as a possible mechanism of electron pairing in 
superconductors (see, i.e., Ref 8 and references therein). 
Nevertheless, to the best of our knowledge this kind of 
acoustic excitations has not been shown to exist experi- 
mentally. Very recently detailed ab initio calcul ation s of 
the dynamical dielectric response of bulk MgBsP^ and 
PcP21 predicted the existence of acoustic plasmons. In 
particular, in the case of Pd this acoustic mode presents 
a 3D character as it is found in all three high-symmetry 
directions!^ 

The main aim of the present study is twofold. First, 
the complexity of the low-energy electronic response of 
bulk Pb is shown, together with the effect of the different 
physical ingredients. Second, the results on the acoustic 
modes are presented and analyzed in detail. Low-energy 
collective excitations show a strong anisotropic character. 
The calculations of the low-energy collective excitations 
have been done using the first-principles pseudopotential 
approach within the time-dependent density functional 
theory (TDDFT)PEl 

The rest of the paper is organized as follows: In Sec. 
II the details of the ab initio calculations of the dynam- 
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FIG. 1: (Color online) Calculated band structure of bulk lead, 
with (solid lines) and without (dashed lines) spin-orbit cou- 
pling in the Hamiltonian. The horizontal dashed line repre- 
sents the Fermi level. 



ical dielectric response of 3D solids are presented. In 
Sec. Ill the general results for the low-energy regime are 
presented together with comparison with experimental 
optical data, while in Sec. IV. the acoustic-like excita- 
tions are analyzed in detail. Finally, the main conclusions 
are presented in Sec. V. Unless otherwise stated, atomic 
units are used throughout, i.e., e 2 = h — m e = 1. 



dissipation theorem to the dielectric function e(r, t',lj). 
For a periodic solid, 
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Im[e G 1 =0)G/=0 (q,o;)] , (1) 



where Q is the normalization volume, G's are reciprocal 
space vectors, and Im[£ G ;L G / =0 (q, <*;)] is the so-called 
energy-loss function, whose Fourier coefficients are re- 
lated to those of the density-response function x(r, t',lj) 
through 

^G' fa' W ) = S G,G> + ^G (q)XG,G' (q, ^) , (2) 



where ^g(q) 



47T 



is the Fourier transform of the 



|q+G|2 
Coulomb interaction. 

The key quantity in the evaluation of Eq. ([2| is the 
density-response function, which in the framework of 
TDDFTE^ll^ satisfies the matrix equation 



xg,c (q, w) = xg,g' (q> + ^2 x g,g- (q, ^) 

G" G"' 

X [v G » (q)*G",G'" + #G",G'" (q)]XG"',G' (q, u), (3) 



II. CALCULATION METHOD 

A crucial quantity in the study of the dynamical dielec- 
tric response of solids is the dynamical structure factor 
5(q, uS) as, within the first Born approximation, the in- 
elastic scattering cross section of X-rays and electrons is 
proportional to it. *S(q, uj) is related by the fluctuation- 



where Xg G' (q> °°) are the Fourier coefficients of the 
density-response function of noninter acting Kohn-Sham 
electrons, i^ccfq) stands for the Fourier components 
of the exchange-correlation kernel, whose exact form is 
unknown. Thus, approximations have to be used to de- 
scribe i^G,G'(q)- 

Instead of directly calculating x GG ,(q, a;), we calcu- 
late first the spectral function matrix 5q g , (q, cj)I 14 * 15 I 



~ BZ occ unocc 

SG,G'(q»") = n EE E (V'nk|e- i(q+G) - r |^ k+q )(V'n'k +q |e i(q+G ' ) - r |^„ k )«5( £nk - e n , k+cl +co) . (4) 

k n n' 

I 



In Eq. Cj), the factor 2 accounts for the spin degeneracy, 
n and n cire band indexes, k is in the first Brillouin zone 
(BZ) and £ n k and i/j n k are Bloch eigenvalues and eigen- 
functions, respectively, of the Kohn-Sham Hamiltonian. 
From the knowlegde of 5 G G ,(q,u;), the imaginary part 
of x G G , (q, uj) is easily evaluated using the expression 

Sc.G'Cq^) = --^n(cj)Im[xG, G '(q^)] > ( 5 ) 

where sgn(uj) = 1 (— 1) for lj > (u < 0). The real 
part of Xg g' (q> u ) 1S obtained from the corresponding 
imaginary part using the Hilbert transformation. 

For G^G', the dielectric response couples the contri- 



bution of q + G 7^ q + G' that appears as a consequence 
of the existence of variations of the electronic density in 
real solids, see Eq. Q. These couplings represent the 
so-called crystalline local-field effects (LFE). If LFE are 
neglected, the energy-loss function is simply given by 

Im [eG=G'=o(q>^)] = 

k G =G'=o(q,^)l 2 

= Im [eG=G / =o(q ? ^)] ^ 

[e G=G / =0 (q, a;)] 2 + Im [e G=G / =0 (q, uj)} 2 

For the density functional theory (DFT) ground state 
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FIG. 2: (Color online) Energy-loss function surfaces for bulk Pb in arbitrary units, with (a),(d) q along T-X, (b),(e) q along 
T-K, (c),(f) q along T-L. (a)-(c) are the results obtained neglecting LFE, while (d)-(f) include LFE. In all six plots, q runs up 
to the r point of the corresponding second BZ. Vertical black lines mark q=X,K,L. The colour bar applies for all the plots 
except the insets. For each high- symmetry direction, q l5 q 2 , q' x and q(> mark the fixed values of the momentum transfer for 
which the energy- loss function is shown in Figs. 4][6 uj in Ry and q=|q| in Bohr -1 in all plots. 



calculations, the electron-ion interaction is represented 
by norm-conserving non-local pseudopotentialsp^ and 
the local density approximation (LDA) is chosen for the 
exchange and correlation potential, with the use of the 
Perdew-Zungei^ parametrization of the XC energy of 
Ceperley and Alder Well-converged results have been 
found with a kinetic energy cut-off of 14 Rydberg (Ry), 
including ~ 150 plane- waves in the expansion of the 
Bloch states, for the face-centered cubic bulk Pb with 
the experimental lattice parameter a = 4.95 A. 

Two different sets of calculations were carried out in 
evaluating Eq. Q. First, ^QQ/fq, uj) was calculated in 
the range < uj < 80 eV with a step of Auj = 0.04 
eV, the band indexes in Eq. Q running up to n = 
25. A Monkhorst Pack 144x144x144 grid of k vectors 
was used as the Brillouin zone sampling with w 32000 
vectors in the irreducible BZ (IBZ). The delta function 
was represented by a Gaussian of width 0.05 eV. Second, 
Eq. ([§ was evaluated for very-low energies, in the range 
< uj < 4 eV, with a step Auj — 1 meV, summing up 
to 12 bands. In this second set of calculations a finer 
432x432x432 grid was used, with w 850000 k vectors 
in the IBZ, and the width of the Gaussian replacing the 
delta function was 2 meV. 

In Fig. [I] the calculated band structure along high- 
symmetry lines of the first BZ of bulk fee lead is shown, 
including (solid lines) and excluding (dashed lines) the 
spin-orbit interaction in the Hamiltonian of the sys- 



tem. As the fee is centrosymmetric, due to the Kramers 
degeneracy^ each band is double-degenerated in both 
cases. The calculated band structure is in good agree- 
ment with other theoretical references (see, i.e., Refl20|) 
and with the experimental one^ when the SO term is 
taken into account. As can be seen in Fig. [TJ the 
inclusion of the SO interaction has sizeable effects on 
the bands crossing the Fermi level (of p-orbital charac- 
ter), mainly around the high-symmetry reciprocal-space 
points. Thus, SO effects on the band structure of bulk 
Pb are expected to show up on its low-energy dielectric 
properties, see Fig. [3j 

In the present work, in the band structure calculations 
the Hamiltonian including the SO term was solved fully 
self consistently. Then, in calculating *S f Q G /(q, uj) SO 
coupling enters Eq. Q only through the energy spec- 
trum via the 5 function. The coupling matrices (the 
brackets in Eq. Q) are kept the same (the ones with the 
wavefunctions calculated at the scalar-relativistic level) 
when using the scalar-relativistic energy spectrum or the 
one with the SO splittings included. Since SO-split bulk 
band structure of Pb, that possesses inversion symmetry 
is reproduced in the first order perturbation theorypl 
we only expect a weak modification of scalar-relativistic 
wave function by SO coupling. Therefore, we do not ex- 
pect that the replacement of the spinor wave function by 
a scalar one should significantly modify the evaluated SO 
results. 
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FIG. 3: (Color online) Comparison of calculated and exper- 
imental optical data. Theoretical results correspond to the 
smallest q G T-X, with q=|q| ~ 0.009 Bohr" 1 . LFE are 
included and the RPA is used. Experimental data at both 
temperatures from Ref. 1241 



III. GENERAL RESULTS 

The simplest estimation of the plasma frequency for a 
bulk metal is given by the FEG model as uj p — ^Ann with 
n being the electronic density, which in terms of the den- 
sity parameter reads uo p For lead, using: the ex- 
perimental value of r e s xp = 2.298, one gets 00p EG = 13.53 
eV. This theoretical value is in good agreement with the 
one obtained in electron energy- loss experiments, see ref- 
erences in Ref. [23] Thus, the energy transfer range uo ^8 
eV is well below the main bulk plasmon energy in bulk 
lead. Note also that calculating the dielectric response 
up to 80 eV the results are well converged with respect 
to the finite energy range used in the Hilbert transform. 

In Figure [2] the energv-loss function surface 
Im [eQ^ G / =0 (q, calculated using the TDLDA ap- 
proximation (see below) and including the SO-induced 
energy splittings is presented as a function of both the 
momentum transfer q and energy transfer cj, with q 
belonging to three different high-symmetry directions, 
both neglecting and including LFE. Note the strong 
anisotropic character of the dielectric screening. In all 
three directions of q several peaks raised by interband 
and intraband transitions disperse. However, there are 
important differences. First, only for q G T-L there is 
one well-defined interband peak, see Fig. [4] Second, the 
number of interband peaks is higher for q G T-L than in 
the rest of the reciprocal space directions. Nevertheless, 
the main anisotropy effect is the distinct topology of 
the energy- loss function for q G T-K. The dispersion of 
the main interband peak is not flattened for q in the 
BZ boundary. And interesentingly, for q vectors in the 
second BZ the dielectric response it is structureless, with 
no feature raising on top of a roughly flat background. 

Comparison of the calculated energy-loss function with 
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FIG. 4: (Color online) Fixed-q Im[s~^ G / =0 (q, uj)] curves, 
q along T-X in (a) and (b), along T-K in (c) and (d), and 
along T-L in (e) and (f). 1BZ and 2BZ stand for first BZ and 
second BZ results. Black lines correspond to q close to the 
first or second BZ center [qi and q[ in Fig. [5](a)-(c)], whereas 
light grey (green in online version) lines correspond to q close 
to the BZ border [q2 and q' 2 in Fig. [2] (a)-(c)]. Solid lines 
are results obtained with inclusion of LFE and SO-induced 
splittings at the TDLDA level (see text). Dashed curves are 
the corresponding results neglecting the LFE. 
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FIG. 5: (Color online) Same as Fig. [4] but with the dashed 
lines showing the energy-loss function calculated with inclu- 
sion of LFE and SO-induced splittings at the RPA level. 



optical experimental data from Ref. [24] obtained at 140 
K and room temperature is shown in Fig. [3] The cal- 
culated curves correspond to the smallest q used along 
T-X, with q=|q| ~ 0.009 Bohr -1 , allowing the compar- 
ison with optical experiments. As can be seen from the 
figure, including the SO-induced splittings in the band 
structure improves the overall agreement with the exper- 
iments, specially with the measurements carried at 140 
K. The deviations of the energy transfer positions of the 
different features with respect to the experimental ones 
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decreases when including the SO-induced splittings in the 
calculation. However, the strong increase of the intensity 
of the interband peak at ~1.7 eV when including the SO 
splitting seems to worsen the agreement with the exper- 
imental measurements. Note that in the experimental 
data of Ref. [24] the two first peaks measured at T=140 
K are washed out at room temperature. Thus, one ex- 
pects the first interband peak to increase its intensity 
when measured at lower temperatures. This is in agree- 
ment with the fact that Pb presents a strong electron- 
phonon coupling, which should broaden the one-electron 
energy levels with increasing temperature. Then, taken 
into account that the methodology used in the present 
work applies for T=0 K, the conclusion is twofold: first, 
a satisfactory agreement with experiments is found which 
supports the approach used in the present work. Second, 
the comparison with the experimental optical data shows 
clear SO effects. 

Below, a one-by-one systematic analysis of the effect 
of the main physical ingredients on the low-energy elec- 
tronic collective excitations in bulk Pb is presented. 



A. Local-field effects 

As shown in Fig. |2j LFE affect differently the cal- 
culated dielectric response of bulk Pb depending on the 
direction of q. As a common trend their effect is stronger 
for q vectors in the second BZ. As can be seen from Fig. 
[4j neglecting LFE decreases the intensity of the energy- 
loss function for all three directions of q, except for small 
q values (compare solid and dashed black lines in Fig. [4]). 
More precisely, the two main effects of neglecting the LFE 
are washing out the strong structureless background for 
q in the second BZ in T-K, and decreasing the intensity 
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FIG. 6: (Color online) Same as Fig. [4] but with the 
dashed lines showing the energy-loss function calculated at 
the TDLDA level with inclusion of LFE and using the scalar- 
relativistic energy spectrum, without SO-induced band struc- 
ture splittings. 



of the main interband peak in T-L as q increases. The 
latter peak is shifted down in energy as much as 0.4 eV 
when neglecting the LFE. 

B. XC kernel 

When calculating the density-response function 
x(r, r',uj) of the interacting electrons from the knowl- 
edge of its non-interacting counterpart, short-range ex- 
change and correlation (XC) effects are accounted for 
by the so-called exchange-correlation kernel, K xc . The 
exact expression of K xc is unknown, and one has to 
choose an approximation. Here, two different approaches 
are used for K xc : the random-phase approximation 
(RPA) and the time-dependent local-density approxima- 
tion (TDLDA) 25 , an adiabatic extension of the local- 
density approximation. In the RPA one simply neglects 
the short-range exchange and correlation effects, thus, 

In Fig. |5j the energy-loss function for several values of 
q belonging to the three different high-symmetry direc- 
tions is plotted, showing results obtained with both RPA 
and TDLDA. As can be seen, the effect of the TDLDA 
with respect to the RPA is the increase of the intensity 
of the calculated Im [^ =G f =0 (^ oj)] , but without chang- 
ing its shape. Only the main interband peak in T-L is 
slightly shifted down by ~0. 1-0.2 eV, see Fig{5](c). 

C. SO-induced splittings effects 

SO-induced energy band splittings also affect the di- 
electric response of bulk Pb in an anisotropic way. As 
shown in Fig. [6j in all three high-symmetry directions 
SO splittings lower the intensity of the broad feature at 
lj ^6 eV. However, including the SO splittings in the cal- 
culation of the energy-loss function has sizeable effects for 
q G T-X and T-L, but it is not the case for q G T-K. In 
T-X, SO-induced band structure splittings broaden the 
scalar-relativistic interband plasmon around the first BZ 
border, shifting it to higher energies by as much as 0.7 
eV. For q G T-L, the main interband peak is shifted up 
in energy by ~0.4 eV upon inclusion of the SO splittings 
in the calculations. 



IV. ACOUSTIC EXCITATIONS 

In the insets of Fig. [2] (d)-(f) the energv-loss function 
for (q— ^ 0, lj —> 0) is zoomed in. In all three direc- 
tions of q linear dispersing features which present van- 
ishing energy transfer at vanishing momentum transfer 
appear. To study the exact character of theses acoustic- 
like modes, highly detailed calculations are needed in or- 
der to resolve the thinest details. Thus, results obtained 
using the second set of calculations (see Sec. II) are an- 
alyzed in the present section. 
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FIG. 7: (Color online) Imaginary and real part of the dielec- 
tric function of bulk Pb, as well as the corresponding energy- 
loss function. All curves correspond to the same q=|q| ~0.027 
Bohr -1 , with q G T-X (solid lines), q G T-K (dashed lines) 
and q G T-L (dashed-dotted lines). The dashed-dotted-dotted 
curves correspond to the Lindhard dielectric function for bulk 
Pb r e s xp = 2.298. The vertical dashed lines (labeled A,B,C) 
mark values of uj where Re[e]=0 for q along T-K (see text). 
All in arbitrary units. 



As can be seen from the insets in Fig. [2j the acous- 
tic modes in all three high-symmetry directions disperse 
keeping its character up to q~0.14 Bohr -1 and uj ~2 eV 
(~0.14 Ryd in the units of uj in Fig. [2|, when they get 
broaden out or merge with interband peaks. All modes 
also present a similar velocity of dispersion. 

In Fig. [7| the calculated imaginary and real part of 
the dielectric function are presented, together with the 
corresponding energy-loss function with q along all three 
different high- symmetry directions. The results were ob- 
tained including LFE, SO splittings and the TDLDA ker- 
nel in the calculations. The Lindhard dielectric function 
for r e s xp = 2.298 and the energy- loss function derived 
from it are also shown (dotted lines). All curves corre- 
spond to q=|q| ^0.027 Bohr -1 . Comparison with the 
Lindhard screening (which stands for the dielectric re- 
sponse of a FEG) shows important band structure effects 
on the calculated e(q, cj), which are responsible for the 
appearance of the acoustic modes. 

As previously proven, in the range (q— >• 0, uj — » 0) 
neglecting LFE has little impact on the results (compare 
the black dashed and solid lines in Fig. [3]). This allows 
for the analysis of the results shown in Fig. [7] to be 
based on the expression for the energy-loss function with 
no LFE, Eq. [6] Thus, if LFE can be neglected, peaks in 
the energy-loss function correspond to minima on both 
the real and imaginary parts of e coinciding at the same 
(q, uj) values. This is in constrat to the definition of 
plasmon in classical electrodynamics^, which states that 
such collective excitations occur whenever Re[e]=0. 

In Fig. [7] the dashed vertical lines mark the energy 



transfer values at which Re[e]=0 for q G T-K. As for 
the Lindhard screening case, there is one first Re[e]=0 
crossing labeled by A, as in all the rest of curves. But 
only for q along T-K Re[e] crosses again zero, at B and C. 
At B, Im[e] has a local minimum, whereas at C it presents 
a local maximum. As a consequence (see Eq. |6|, at B a 
peak rises, which is not the case at C. 

For q G T-L, Re[e] crosses zero only once, but it is close 
to vanish around ~0.35 eV, ~0.4-0.45 eV and ^0.55-0.6 
eV. In all three cases, there is a local minimum in the 
imaginary part of the dielectric function. As a result, 
three peaks raise in the corresponding energy-loss func- 
tion (see the green line in the lower panel of Fig. [7]). 
However, from the point of view of classical electrody- 
namics, the character of the acoustics modes for q along 
T-K and T-L are different as only in the former case the 
real part of the dielectric function vanishes several times 
in the very- low energy transfer range. This might reflect 
the fact that two and not only one band crosses the Fermi 
level in T-K, and thus be related to the acoustic plasmons 
(see the Introduction). 

The vp of the two bands crossing Ep along T-K 
are 14.6 eV-Bohr and 13.3 eV-Bohr. Comparing them 
with the group velocity of the acoustic plasmon in T- 
K, v = 12.7(13.2) eV-Bohr including (excluding) SO cou- 
pling, one sees that v is comparable to the Fermi velocity 
of the bands. A detailed analysis of Im[e] revealed that 
the acoustic plasmon is exclusively the result of the intra- 
band contributions of the two bands crossing Ep along 
T-K, as expected for Im[e(q— » 0, uj — » 0)]. 

The difference between the modes along T-K and along 
the rest of high-symmetry directions may not only be a 
matter of classification. Taking into account that the 
methodology used in the present work applies only at 
T=0 K, following the discussion of the comparison with 
experimental data (see Sec. Ill), at finite temperature 
one would expect the acoustic plasmon along T-K to be 
more suited than the modes along T-L for being resolved 
in electron energy-loss experiments. If the modes with 
q G T-L broaden out, they would easily overlap with each 
other resulting in a highly broad feature. On the other 
hand, the well-defined acoustic plasmon in T-K presents 
a better background-peak intensity ratio and a smaller 
linewidth. 



V. CONCLUSIONS 

We have presented first-principles calculations of the 
low-energy (uj ^8 eV) electronic collective excitations 
in bulk Pb and studied in detail the effect of the main 
physical ingredients involved, as well as the existence and 
character of acoustic-like modes. Good agreement with 
available optical experimental data is interpreted as an 
evidence of strong SO effects. 

In general, strong anisotropic effects are found, results 
showing a distinct topology of Im[^Q^ G /= (q, oj)} for 
q G T-K. LFE and SO-induced splittings have sizeable 
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effects on the dieletric screening of bulk Pb, showing an 
anisotropic behavior. For q vectors in the second BZ, 
LFE on the energy-loss function are remarkable. Inclu- 
sion of exchange-correlation effects through the TDLDA 
kernel increases the intensity of the energy-loss function 
in the studied range, but without affecting its shape. 

Very-low energy modes with acoustic-like dispersions 
are found in all three studied high-symmetry directions 
of q and are shown to be a consequence of band structure 
effects. The character of these acoustic modes depends 
on the direction of q. The possibility of experimental 
check of the existence of these acoustic modes by elec- 
tron energy-loss measurements seems feasible as these 
modes keep their character up to uj ~2 eV. And as a 
result of a detailed analysis, the acoustic plasmon along 
T-K is suggested as the most suitable for its experimental 
detection. 

As an overall conclusion, the present study has shown 
the complexity of the low-energy dielectric response of 
bulk Pb, even at the (q— >• 0, uo — >• 0) limit through 
the anisotropic nature of the acoustic modes. Thus, 



in theoretical studies of more complex properties (i.e., 
two-band superconductivit}^), care should be taken 
in using models for the dielectric properties of fee bulk 
lead, the present results suggesting the need of an ab 
initio description 27 whenever possible. However, in an 
ab initio study on the inelastic lifetime of hot quasipar- 
ticles in bulk leacP^ -which presented agreement with 
experiment j 29 l 3Q l . no remarkable dependence on their 
momentum was found, with a small scattering of the 
calculated linewidths as a function of the quasiparticle 
energy (see Fig. 3 in Ref. |28|) . 
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